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Abstract On the basis of the study of inclusive hadronic events, two methods are 
adopted to determine the number of produced ip{2S) events collected by BES in 2001—2002 
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run, which is 14.0 x 10 6 with the uncertainty of 4%. 
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1 Introduction 

The BEijing Spectrometer (BES) is a general purpose solenoidal detectorE running at Beijing 
Electron Positron Collider (BEPC). The beam energy of BEPC is in the range from 1.5 GeV to 
2.8 GeV with a design luminosity of 1.7 x 10 cm~ 2 s at 5.6 GeV center of mass energy. The 
main physics goal is to study the charm and r physics. During 2001-2002 years' running, about 
14 million ip(2S) online hadronic events have been collected.]. On the basis of this large data 
sample, many physics analyses could be performed with an unprecedented precision. 

The determination of the offline total number of ift(2S) event, N^2g^, is a foundational work 
in physics analysis, and in turn is the foundation of the further analysis study. In tp(2S) physics 
analysis, the calculation of the absolute branching ratio depends on A^J^ , whose error will be 
directly accounted into the error of the branching ratio of any being studied channel. Therefore, 
it is essential to work out the A^J) accurately and reliably. In principle, any decay channel 
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totally 3000 runs (RUN20050-RUN23085) are taken, among which only those with I qua ut y = 2, 3 are 
used to determine the total number of the ip(2S) events. Here I qua iity denotes the run quality and value 
2 and 3 indicate that the run quality is fairly reliable. 



with known branching ratio could be used to evaluate the total number of ip(2S): 

at tot _ iV / 

N *w ~ TPs] ' 

where Nj bs is observed number of final state / for a certain decay channel, ej and 13 f are the 
corresponding efficiency and the branching ratio. It is obvious that the larger the branching 
ratio and the smaller the corresponding error, the more reliable the total number is. On such an 
extent, the inclusive hadron final state is a favorable process for the total number determination. 
The only disadvantage here lies in the difficulty to eliminate all kinds of backgrounds throughly. 
Therefore the meticulous studies have been made for the hadron event selection. 

In the following sections, the hadron event selection is discussed firstly, then two methods 
are utilized to determine the total number of ip(2S) and the uncertainties from various sources 
are studied. At last, the final result is given. 

For clearness and convenience, some notations which are to be used afterwards, are listed in 
the TableE 

Table 1: Notations 



Symbol 


Meaning 


Superscript 


Meaning 


Subscript 


Meaning 


N 


Observed Number 


T 


Total 


h 


hadron final state 


m 


Selected Number 


P 


Peak region 


e 


e + e~ final state 


n 


"Pure" Hadron Number 


R 


Resonance 






e 


Efficiency 


C 


Continuum 






a 


Cross section 










L 


Integrated Luminosity 











In addition, there are two elementary relations among the five quantities N, n, e, a, and L, 
that is 

N 

L= — or N = L-a , (I) 
a 

Tl 

N = — or n = N ■ e . (2) 
e 

The symbol with a tilde on it (e.g. h, N, etc.), denotes the events obtained at the continuum 
region (Ebeam = 3.665 GeV), while others denote the events obtained at the resonance region 
(E beam = 3.686 GeV). 

There are also two frequently used equalities, the first one for variables of the same process 
at different energy points: 

N/a L n/(e-a) n L ■ a 

mT = T = ~T( \ , or - = ( for e = e ) ; (3) 

iV/er L n/(e-a) n L ■ a 



the second one for variables of different process at the same energy points: 



n J /e J 



N J 



W 
L-a J 



a 



a-' 



or 



n 



rv' 



(4) 



2 Hadron Event Selection 

For the hadron event selection, the detail information could be found in Refs. [2] and [3]. There 
is no particular event topology to require; instead cuts are made to reject major backgrounds: 
cosmic rays, beam associated backgrounds, two-photon process (7*7*), mis-identified "hadron" 
event from QED processes of e + e~ — > l + l~ , I = e,fi, r, and e + e~ — ► 77 followed by 7 conversion, 
and so forth. Most of these kinds of event have salient topology and could be eliminated by 
proper criteria. Events with at least two well reconstructed charged tracks within | cos#| < 0.8 
are selected (that is N goo( i > 2). The total energy deposited by an event in the BSC (E sum ) 
is required to be larger than 0.36 Eb eam , in order to suppress the contamination from two- 
photon processes and beam associated backgrounds. Events with all tracks pointing to the same 
hemisphere in at least one of axial directions (x or y or z direction) are removed to suppress 
beam associated backgrounds. (This requirement could be expressed quantitatively as I ss i > 1, 
where I ss i is called the squared spatial distribution index.) For two-prong events, two additional 
cuts are applied to eliminate possible lepton pair backgrounds. The number of photons must be 
greater than one (that is N rea i 7 > 2), and the acollinearity between two charged tracks, aAcol, 
must be greater than 10 degrees. 

After the event selection, the fitting of double Gaussian plus a polynomial is applied to 
eliminate the remained background from beam associated backgrounds, see Fig. Q] and [2j 
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Figure 1: Vertex-Fit distribution Figure 2: Vertex-Fit distribution 
[E cm = 3.686 GeV). [E cm = 3.665 GeV). 



The global error analysis approach is adopted to obtain the uncertainties of selection cuts [3] , 
as listed in Table [2) 



Table 2: Error of hadron event selection. 



Requirement 


Error 


Ng 00 d > 2 


2.765 % 




2.544% 


Issi ^ 1 


0.173 % 


Nrea h > 2 (for N good = 2 ) 


0.042 % 


ctAcoi > 10° (for N good = 2 ) 


0.157 % 


Sum 


3.765 % 



In fact, there are many processes that could lead to hadron final state at ip(2S) resonance 



region, they could be divided into seven categories: 

e + e~ — > hadron (CH) , (5) 

e + e" -> ^(25) -> /mdron (RH) , (6) 

e + e~ — ► J/V> — > hadron (J/ip-H) , (7) 

e + e~ -> r + r~ Ziaciron (r-CH) , (8) 

e + e~ -» if>(2S) -> r + r~ -> Ziadron (r-RH) , (9) 

e+e" -» hadron* (CH*) , (10) 

e+e- -> ^(25) -» Ziadron* (RH*) , (11) 



where C represents the continuum process, R the resonance process, H hadron event, and H* 
indicates the event which survives all aforementioned hadron selection cuts and is left in hadron 
sample from processes (|10|) and pip . Among above seven categories, only the hadron event 
from the first two processes are "pure" hadron event at ^{"IS) peak region while others should 
be treated as backgrounds. Since hadron event from different process has almost the same event 
topology, the theoretical estimation method is used to evaluate the contamination of such kinds 
of backgrounds. According to the analysis in Ref. [3], two factors are introduced to subtract 
the hadronic background. If the pure hadron number is denoted as n and the selected hadron 
number denoted as m, then it could be obtained 

1cm^ ad = n% ad , (12) 
iR^had = n had , (13) 

where 

ic = 1 - f JH > - fc = 0.855 , 
1R = l-f R = 0.999 , 

with 

f had had 

JJ/i> ~ c . f c ' 

had had 



I 



fc 



E 



C JC 



j9 ■< 

had had 



C 



Ir 



E 



B 



I 



Here all efficiencies are obtained from Monte Carlo simulation [3], and symbols r, e,/x,7 and 7* 
denote final states r + r~, e + e _ , fi + fi~, 77 and 7*7*, respectively. For the continuum process, 
the production cross section could be obtained from the corresponding Monte Carlo gener- 
ator. For cr? ad and <7^\ they could be calculated by theoretical formulae with corresponding 
resonance parameters obtained from the scan experiment data which is about 15 nb and 

1 nb, respectively. Since Monte Carlo generator does not give the cross section for resonance 
process, the ratio of branching fractions is used in the calculation of factor /r [SJ. It should 
be pointed out that because the variation of is fairly smooth at ip(2S) region, refer to 
Fig. [3j the contamination from J ftp decay is suitable to be treated as the continuum- like hadron 
background. 



10 



10 



lO 




3.1 



3.2 



3.3 



3.4 



3.5 



3.6 



3.7 



Figure 3: Cross section curve of J/ip resonance. 



3 Determination of Total Number 
3.1 Principle 

The branching ratio of hadron final state, denoted as Bh, can be acquired from PDG list [5] 
or from BES scan results [2]. If the number of selected hadron event from ijj^lS) resonance is 
n i>(2S)^had., them 

QT n^2S)^had. , s 

e h °h 



However, the number of selected hadron event at a certain energy point is the combination of 
two parts (refer to Fig. one from resonance process and the other from continuum one, that 
ifl 

nhad = n$(?S)-*kad. + n e+e-^had. > 

4 4 4 (15) 

™h = n h( n h) + ra £ • 
Therefore the key issue here is to distinguish the from the n\. There are two methods, the 
fraction subtraction method and the normalization subtraction method, can be used to figure 
out the number of resonance event. 
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Figure 4: Fractional sketch of hadron number of different process. (The cross-hatched 
area in the figure indicates the data taking region.) 



3.2 Fraction subtraction 

Refer to Fig. HI if the ratio of to could be estimated: 

j p i 

together with Eq. (|15p . the can be calculated as 

As to the factor /, from Eq. (jj]), it is easy to acquire 

n h _ e h ' a h G _ e h • a h P 

„P ~ e R „R ' U1 n h — r r n h ) 
n h e h ' a h e h ' a h 



2 Because of energy spread and shift, the data are actually taken within a region rather than at a 
fixed point (notice the cross-hatched area in Fig. 0|. So "P" instead of "R" is adopted in the following 
formulae in order to denote such an effect. 



a 



so / can be expressed as 

e c ■ rr c 

f = J n ■ ( 17 ) 
e h - G h 

Combing Eqs. (fTf]) and (fTo]) . it can be obtained 

T 

atTOT _ n h /-.ox 

^ i2S) ~ B h -^-(l + f) ■ (18) 
3.3 Normalization subtraction 



The data taken far from the resonance region could be treated as the data of continuurrj§, refer 

l h 



to Fig. U from Eq. ([3]), the number at continuum region n? could be transformed into that at 



resonance by a luminosity normalization factor 



n C h =frn C h , fT = %^ 



where fx is the transformed factor. Combining with the relation = n£ + , the resonance 
number could be worked out 

n h 

In the expression of fx, the luminosity L is usually calculated by the continuum e + e~ event (i.e. 
Bhabha event): 

L 



Similar to Eq. (|16p . at the continuum region, the number of events from the continuum process 
is expressed by the number of total selected events: 

T ~T ~T 



n? = — , n? = — e^r- , h c 



l + /e ' £ ~l + fe ' h ~l + f h 



with 



c e u e 

/e = #T# - (20) 

h = y$ . (2D 
e h ■ a h 

It should be noticed that the relation = has been used 1 in the above calculation. 



3 Strictly speaking, within the scan range, any data are from two processes, resonance and continuum. 
The effect of resonance to continuum could be taken into account by factor / e , / e , fh in Eqs. (|19p . (120]) . 
and ([21]). 

4 The relation = ef\ is exact for the e + e~ final state, whose event selection cuts are energy 
independent; but for the hadron final state, the relation is only an approximation. 



Put all together, the total number can be calculated as 



where 



N ™2S) = A-(1-FT-M eh ) , (22) 



_ a c h -~a c e (1 + h) 1 , , 

* T -tf + fe) '(! + /,) ' ^ 



with factors f e , f e , and fh are given in Eqs. (fTU|) . (|2"0j) and (|2"Tj) . 

3.4 Correction 

As mentioned in section [2j after the hadron event selection, the selected hadron number m 
instead of pure one n is obtained. The relation between m and n is given in Eq. (|12p and (|13p . 
that is 

7c(m£,m£) = (n£,n£) , 
j R (m%,m%) = (n£,n£) . 



Notice that 



then 



that is 



7i? 7C 1R 1C \1R 7c* 



p m 



J_ , X 

7R 7C 

So the formula for the fraction subtraction method now becomes 



N TOT _ fj}h (0 r\ 

N *CX)- Bh .e*-{l + f).&, ' (25) 



'71 



where <5 71 is defined as 



*n = ^f^f 1 • (26) 

For the normalization subtraction method, the corresponding corrected formula could be 
obtained similarly and the final result is 

where 



and 

•Ft = Ft/5 12 
here Ft is just as that defined in Eq. (|23p . 

3.5 Numerical Calculation 

By the virtue of Eqs. (|25|) and (I27p . N^2gs could be worked out. For convenience, all numbers 
relevant to total number calculation are summarized in Table El 

Table 3: Event number, efficiency and cross section. 



5 12 



l + fh 



(29) 



Final state 


hadron 


e + e 


Region 


Resonance Continuum 


Resonance Continuum 


Event number 


ml = 10634586.0 ml = 13937.7 
(Ami = H688.2) (Am[ = 817.2) 


nf = 1190613 hi = 54415 


Efficiency 


Resonance 
process 


el = 0.753 ef « e« 
(Aef = 0.012) 


ef = 0.761 ef = ef 


Continuum 

process 


e c h = 0.716 et « e£ 
(Aef = 0.010) 


e f = 0.708 ef = ef 


Cross 
section 
(nb) 


Resonance 
process 


of = 676.277 of = 0.266 


of = 4.053 of = 0.00158 


Continuum 

process 


of = 15.495 of = 15.673 


of = 80.423 crf = 81.349 


Trigger efficiency fl 


Correction factor 


= 0.99924 


lR = 0.999 7c = 0.855 



f: the trigger efficiency is obtained from Ref. |6J written by Dr. Fu ChcngDong. 



N^2g\ is worked out to be either 14.05 xlO 6 for Fraction method or 14.04xl0 6 for Normal- 
ization method. 



4 Error Analysis 
4.1 Classification 

Formally, by the virtue of Eqs. (|18f) and (|22h . the formula to calculate the total number can be 
written as 



^ = -rV G « ' (30) 



T 



where G{ is a correction factor denned as 




for the fraction method; 



G l = { 1 + / ' 

= (1 — Ft • M e h) , for the normalization method. 

The error of ^VjSJ) comes from the components of Eq. ([30]) . such as n T , eff, B^, and Gj, the 
error of which will be discussed one by one. 

4.2 Uncertainty of selected number 

For the selected numbej§ mT, there are three sources of uncertaintjl^: 

1. Fitting uncertainty 

The uncertainty of fitted number could be obtained from the corresponding error of fitting 
parameters which are used to calculate the number (refer to Table [3|) , that is 

Vfit(rn T h ) = = 0.11% 

m h 

2. Statistic uncertainty 

According to statistic principle, 

Vstaijnl) = = 0.03% 

y m h 

3. Selection uncertainty 

According to the study of section 1, the selection uncertainty reflects the inconsistency 
between data (m^) and Monte Carlo (e^ ), so the uncertainty of is also included in this 
term which is 

^ Of) = 3 ' 77% " 

4. Effect due to beam energy fluctuation 

As it is mentioned in footnote [5J the data are actually taken within an enegy range (refer 
to the sketch description in Fig. H]), so the ratio between cr^/cr^ will vary with the actual 
beam energy which may be different for different beam-injection. Usually, each beam- 
injection includes 3—5 runs. As an estimation, all runs are grouped with every 3, 4 or 5 
runs, then the ratios of selected hadron number (n^) to that of e + e~ number (n£) are 
worked out, which denoted as r^i) with i indicating the grouped run number. Fig. [5] 
shows the r^e distribution for different grouped-run number. The maximum value of rh e (i) 
corresponds to the peak cross section^. Taking the experimental statistic fluctuation into 



5 Hereafter the selected number m instead of pure number n is used in the error analysis, and the 
uncertainty for such substitution is rather small and is to be discussed afterwards. 
6 Hereafter the symbol v denotes relative error. 
7 Notice 

rih , T n e 
Oh = t , and L = , 



i n 



account, the value 10 is adopted as the position of the peak cross section. The uncertainty 
from the beam energy fluctuation effect is estimated as follows 



'ex~p 



</(! + /) 
5>£(i)/[l + /(•)] ' 



with 



10 



For different grouped-run cases, u exp is almost the same, which is 0.23 %. 
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Figure 5: The distributions for different grouped-run number. 



4.3 Uncertainty of branching ratio Bh 

There are two values, one from PDG2002 and the other from ^{IS) scan experiment: 

£/ l (PDG2002) = (98.10 ±0.30))% , 
B h (^(2S) scan) = (97.79 ± 0.15)% . 

In the total number calculation the Bh value form PDG2002 is used. However as a conser- 
vative estimation, the difference between above two values is used as the uncertainty of Bh, that 
is 

u{B h ) = 0.32% . 

then 

CT e • e e n h n h 
a h = • — , or a h oc — . 

e/i n e n e 



4.4 Uncertainty of correction factor Gi 

For G\, the uncertainty due to different / could be calculated as follows 



Ki2 T s) u) 



1 + / 1 + /' 



1 



1 + / 



/'-/ 



1 + f 



According to the definition of /, Eq. (|17|) . the uncetainty of / mainly comes from the statistics 
of e£ and , both of which are 1/V50000. Notice / = 0.02179 is a small quantity and the 
maximum difference between / and /' is also small, so 

i/(Gi)«|l-/|-A«0. 



where A = \f - f'\ max = 2/^50000. 

For Gi-, the uncertainty due to different G 2 could be calculated as follows 

*Jg§(G3)-JVjgSS(Gi) 



u(G 2 



F T ■ M eh - F' T ■ M' eh 



I -F T - M, 



e/i 



^gJ)(G 2 ) 

Notice Ft{= 0.932) approximates to one, so the difference of total number calculated with 
Ft = 1 and with Ft 7^ 1, is treated as the error from factor Ft, that is 

M eh ■ (F T - 1) 



v(F T ) 



0.20% 



I -F T - M eh 

Notice M e h consists of two pairs of ratio, so the systematic error of numerator and denominator 
will cancel automatically, only the statistic error is left, which is 



AM 



1 



1 



1 



1 



eh 



~p I ~ i~p I rp I 1 i 

nj nj nl n L h 



0. 



so the error due to different M e u could be calculated as 



Then the uncertainty for G2 is 

u(G 2 



F T ■ (M eh - M' e 



eh) 



I -F T - M eh 



< 



F T ■ AM, 



eh 



I -F T - M eh 



0.91% 



v 2 (F T ) + v 2 {M eh ) = 0.94% 



4.5 Other uncertainty 

The other effects which could lead to the uncertainty include the correction factor 7c, jr, trigger 
efficiency and so forth. All these uncertainties are regarded as less than 0.5%. In addition, Monte 
Carlo sample is used as input to test the bias of our method. The error from such bias is about 
0.6%. 



4.6 Summary 

Put all things together, the synthetic uncertainties are summarized in Table HI 

1 o 



Table 4: Error summary. 



Source 


Uncertainty 


T 


Fitting 


0.11 % 




Statistic 


0.03 % 




Selection 


3.77 % 




Ebeam fluctuation 


0.23 % 




0.32 % 


Gi 


Gi 


0.88 % 




G 2 


0.94 % 


Other 


(0.500.6) % 


Total 


Method 1 


3.97 % 




Method 2 


3.98 % 



Note: Method 1 for Fraction method; Method 2 for Normalizaion method. 



5 Result and Discussion 

The final total number of ip(2S) event with corresponding error is 

TOT _ ( 14.05 x (1 ± 3.97%) x 10 6 , (Fraction Method) ; 
V(2S) - | UM x q ± 3.98%) x 10 6 ( (Normalization Method) . 

Notice Eqs. (JTHJ) and (|22jl . two methods, the fraction method and the normalization method, 
correlated closely with each other and the difference between two methods for the central value is 
actually less than one per mille. Furthermore, the difference of the uncertainty for two methods 
is also fairly small. Therefore, the central value of the offline total number of ip{2S) event could 
be regarded as 14.0 x 10 6 , and the uncertainty could conservatively be estimated as 4%. 



Thanks are due to Prof. J. C. Chen and Prof. D. H. Zhang for their help about Monte Carlo 
simulation. 
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